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The quantum spin-1/2 two-leg ladder with an anisotropic XYZ Heisenberg intra-rung interaction and Ising 
inter-rung interactions is treated by means of a rigorous approach based on the unitary transformation. The 
particular case of the considered model with X-X intra-rung interaction resembles a quantum compass ladder 
with additional frustrating diagonal Ising interactions. Using an appropriately chosen unitary transformation, 
the model under investigation may be reduced to a transverse Ising chain with composite spins, and one may 
subsequently find the ground state quite rigorously. We obtain a ground-state phase diagram and analyze the 
interplay of the competition between several factors: the XYZ anisotropy in the Heisenberg intra-rung cou- 
pling, the Ising interaction along the legs, and the frustrating diagonal Ising interaction. The investigated model 
shows extraordinarily diverse ground-state phase diagrams includingseveral unusual quantum ordered phases, 
two different disordered quantum paramagnetic phases, as well as discontinuous or continuous quantum phase 
transitions between those phases. 

Key words: quantum spin ladder, exact results, ground state, quantum phase transitions 
PACS: 05.30.Rf, 75.10.Jm 



1. Introduction 

Quantum spin ladders with frustrated interactions are intensively studied during the last few decades, 
since they exhibit a rather complex ground-state behaviour to be reflected in extraordinarily rich low- 
temperature thermodynamics as well Quite recently, a number of exact solutions have been ob- 
tained for several particular examples of quantum spin-| two-leg ladders (3j-l7|]. The railroad ladder con- 
sidered by Lai and Montrunich [6] has a quite specific configuration of inter-spin interactions, namely, 
the staggering of X-X and Y-Y couplings along the legs is supplemented by the uniform Z-Z coupling 
present along the rungs. An exact solution of this specific quantum spin ladder has been found by adopt- 
ing the method originally developed by Kitaev |8], which proved a striking spin-liquid ground state in this 
quantum spin ladder. On the other hand, the railroad ladder with the uniform Z-Z interaction along the 
legs and the uniform X-X interaction along the rungs has been rigorously solved by Brzezicki and Oles 
I3j-l5|]- To a certain extent, this exactly solved quantum spin ladder can be regarded as an one-dimensional 
analogue of the quantum compass model on a square lattice, which describes the orbital ordering in 
transition-metal compounds (9j. 

In this work, we will examine a more general model of the quantum spin-| two-leg ladder, which 
includes the fully anisotropic XFZ-Heisenberg coupling between spins from the same rung and two dif- 
ferent Ising (Z-Z) interactions between spins from neighbouring rungs considered along the legs and 
across the diagonals, respectively. Nevertheless, it should be mentioned that the investigated quantum 
spin ladder extends our previous exact calculations for the spin-j Heisenberg-Ising ladder with the X X Z 
intra-rung coupling [ 7], whereas the quantum compass ladder investigated in detail by Brzezicki and Oles 
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Ud] also represents a very special limiting case of the investigated model system. The main goal of the 
present paper is to examine the simultaneous effect of two different kinds of frustration: the geomet- 
ric frustration caused by the antiferromagnetic interaction between spins from an elementary triangle 
plaquette and the competition between X-X and Y-Y intra-rung interactions with both Z-Z inter-rung 
interactions. 

The outline of the paper is as follows. In section 2, we define the model and show how to get the 
ground state by a rigorous calculation based on the appropriate unitary transformation. The ground- 
state phase diagram of the spin-i XZ-Ising and XY*-Ising ladders is explored in section 3. Finally, some 
conclusions are drawn in section 4. 



2. Model and solution 

Consider the quantum spin-i Heisenberg-Ising ladder with an anisotropic XYZ intra-rung coupling 
and two different Ising-type couplings, which involve Z-Z spin-spin interactions along the legs and 
across the diagonals of a two-leg ladder (see figure[l}: 

N 

1 = 1 

Here, s" . denote three spatial components a - x,y,z of the spin-i operator, the former subscript j = 1,2 
determines the number of a leg and the latter subscript enumerates the lattice position in a particular leg. 
Apparently, the interaction terms J\,]\,J\ account for the quite anisotropic XFZ-Heisenberg coupling 
between two spins belonging to the same rung, while the interaction terms J2 and h take into consider- 
ation the Ising-type interactions between the nearest-neighbor spins along the legs and across the diago- 
nals of the two-leg ladder. It should be pointed out that the z-component of the total spin S? = s* . + s| , 



Figure 1. (Color online) A schematic representation of the considered Heisenberg-Ising two-leg ladder. 
The intra-rung XFZ-Heisenberg coupling is represented by thick lines, while Z-Z Ising interactions 
along the legs and across the diagonals are displayed by thin solid and broken lines, respectively. 

on ith rung commutes with the Hamiltonian, i.e. [{S*) Z ,H] - 0, and this property allows us to present 
the model in an integrable form. To obtain a ground state of the investigated quantum spin model quite 
rigorously, it is advisable to closely follow the procedure worked out previously for the quantum spin-i 
Heisenberg-Ising ladder with XXZ intra-rung interaction. After the unitary transformation (7|] 
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one may rewrite the Hamiltonian i2.ll into the following pseudospin representation: 



H 
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■ J v s z 



[h (1 + 45^5^.^) + 2/3 (si + S* u+1 )] sl;sl i+l 1 , (2.3) 



which shows the symmetry of the model in a more explicit way. It is quite obvious that only z-components 
of spin operators from the first leg are present in the Hamiltonian 12. 3L which means that s* , are good 
quantum numbers. By contrast, different spatial components of spin operators from the second leg are 
still involved in the Hamiltonian 12.3) and thus, they still represent quantum spins with regard to the 
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presence of two non-commuting parts s* . and s| . of each spin operator. Altogether, the Hamiltonian 
d2.3> can be identified as the Ising chain of composite spins in a transverse field, whereas the values of 
the effective interaction and the effective transverse field locally depend on a particular choice of eigen- 
values of the classical Ising spins .. Following 01, one may also establish the following correspondence 
between new and initial states: 



lIl,il2,/> = l^0 0> = -y=(llMt2,/>-ITl,il2,i», 

v2 

Ih.ihi) = \<f>[ _> = 4= (I U.iU.i) - 1 ll.i \2,i)) , 
V2 



lIwTw> = l0'o> = -F(lluTa,i> + lti,il2,i», 
v2 

lhA<> = \<t>[ +> = 4= h >i + 1 ii.^y)) . 

v2 



(2.4) 



where \cp' 00 ) denotes the singlet state of the zth rung and the other three states \cf>[ ,), \<p[ ) form the 
triplet state of the zth rung. It should be stressed that the square of zth component of the total spin S? 
on the ith rung acquires two different values. It either equals (S?) 2 = for \(f>' > and \(p\ ), or (S?) 2 = 1 
for \<p\ + ). To get the partition function one has to diagonalize the Hamiltonian 12. 3> for all particular 
configurations of , and sum up all contributions in the trace of statistical operator. However, it is quite 
evident from the transformed Hamiltonian d2.3t that the chain decomposes into two independent parts 
whenever two neighboring spins S? . and sf , + , have opposite orientation (i.e., take on different eigen- 
values). In this respect, the composite chain is divided into a set of finite chains of different sizes for 
any chosen configuration of .. Generally, this problem seems to represent a quite formidable task, but 
the ground state of the investigated model can be found quite rigorously using the same arguments as 
given in 01 • Since the ground-state energy of two finite but isolated spin-j Ising chains in a transverse 
field is always higher than the ground-state energy of one unique spin-i Ising chain in a transverse field 
obtained by joining both independent finite chains, the ground state of the model under investigation 
should accordingly correspond only to the uniform configuration of all s^ { . Therefore, one may single 
out only two different uniform configurations with all $? , = | or all 5* . = - \ from which the ground 
state of the Heisenber-Ising ladder can be derived. The effective Hamiltonian H2.3t for the two uniform 
configurations acquires the following form: 



A' 



H± = L 
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The ground state energies per site of both effective Hamiltonians can be exactly calculated using the 
Jordan- Wigner fermionization jiollllll : 



(7l +l/2±/3l). 
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is the complete elliptic integral of the second kind. 

Both Hamiltonians H + and H~ imply a precise mapping correspondence between the spin-| Heisen- 
berg-Ising two-leg ladder and the spin-^ quantum Ising chain in a transverse field, which can be, how- 
ever, characterized by different values of the effective interaction and transverse field. Bearing this 
in mind, one should expect quantum phase transitions of two different types. The first kind of zero- 
temperature phase transitions may correspond to a continuous (second-order) quantum phase transition 
inherent to the transverse Ising chain, which arises for one particular ratio between the effective inter- 
action and transverse field. Beside this, there may also occur discontinuous (first-order) quantum phase 
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transitions whenever a crossing of the lowest-energy levels inherent to both effective Hamiltonians <2.5> 
takes place. In the following two sections, we will illustrate all the aforementioned features of quantum 
phase transitions on ground-state phase diagrams of two particular cases of the model under considera- 
tion. 

3. Ground state of XZ-Ising and XF-Ising ladders 

In this section, we will consider two particular cases of the investigated model <2.H by switching off 
either the y- or z-component of XFZ-Heisenberg coupling (i.e., either ft — or = 0). It is notewor- 
thy that the two aforementioned particular cases represent a direct extension of the quantum compass 
ladder to which the considered model reduces when neglecting the z-component of the Heisenberg 
coupling (/j = 0), one of the two transverse components of the Heisenberg coupling (i.e., either Jf = or 
]\ = 0) and the frustrating Ising interaction across the diagonals (J3 - 0). Furthermore, the problem of 
two-dimensional quantum compass model is quite complex and the exact solution for this model has not 
been found yet. 

Let us first consider all possible phases that may appear in the ground state of the model under in- 
vestigation. Each uniform configuration of 5^ . corresponds to the transverse Ising chain, which has three 
possible ground-state phases. The ground-state phases for all . = — i appear if < . It is worthwhile 
to remark that the ground-state phases belonging to this effective model were thoroughly analyzed in our 
preceding paper [7] and let us, therefore, give here just their definition for the sake of easy reference: 

• Quantum paramagnetic (QPM1) state for |(/* + j\) > I/2-/3I: the Heisenberg-Ising ladder is in 
the gapped disordered state with no spontaneous magnetization. The rung singlet dimers dominate 
on the Heisenberg bonds and in the very special case J2 = h, the ground state factorizes to a set of 
the completely non-correlated rung singlet dimers. 

• Stripe Leg (SL) state for \{J* + j() < h~ h- the Heisenberg-Ising ladder shows a ferromagnetic 
order along the legs and antiferromagnetic order along the rungs, i.e., the magnetizations of two 
chains are oriented opposite to each other ({s* .) = (s* . .) = -(sf .) = + 0). The following 
staggered magnetization can be defined as the relevant order parameter of this phase 

< = ^ !(<<!> -<<«>) = 5 

which undergoes the obvious quantum reduction of magnetization. 

• Neel state for j{Jf + J?) < J2- h'- the nearest-neighbor spins both along the legs and rungs ex- 
hibit predominantly antiferromagnetic ordering. The dependence of staggered magnetization as 
the relevant order parameter is quite analogous to the previous case 

The most fundamental difference between the ground states of the Heisenberg-Ising ladder with the 
XX Z- and XFZ-Heisenberg intra-rung interaction can be found in the phases arising from the uni- 
form configuration with all 5* . = | . While in the former model with the XX Z intra-rung interaction, 
all ground-state phases are classical in their character 01. the emergent ground-state phases of the latter 
model with the more anisotropic XYZ intra-rung coupling display significant quantum features. One 
may indeed identify the following three quantum ground states for a particular case < with all 

c z — I ■ 

~ 2- 

• Quantum paramagnetic (QPM2) state for \\Jf- K \ > I/2 + /3I: the equivalent transverse Ising 
chain H + i2.5) is in the gapped disordered state with no spontaneous magnetization (sf) = and 
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non-zero magnetization < sf > + induced by the effective transverse field. For the initial Heisenberg- 
Ising ladder, one consequently gets the ground state with the prevailing dimer state \<p[_) on the 
rungs. 

• Stripe Rung (SR) state for \\Jf-Ji I < h+h- the equivalent transverse Ising chain exhibits a spon- 
taneous antiferromagnetic ordering with {Sf > = (-1) ' m x # 0. Due to relationships 12.2) , one obtains 
for the Heisenberg-Ising ladder {s^ .) = -<sj = is* ,) = -is* + 0. Thus, the Heisenberg-Ising 
ladder shows an antiferromagnetic order along the legs and ferromagnetic order along the rungs. 
The staggered magnetization as the relevant order parameter in this phase is non-zero and it ex- 
hibits evident quantum reduction of the magnetization given by: 



1 



N 



SR 



1) 



Hh + h) 2 



Ferromagnetic (FM) state for j \Jf-j\\ < - Uz + J?,)'- the effective transverse Ising chain shows a 
ferromagnetic ordering. Due to transformation relationships \2.2\ , all the spins of the ladder have 
the same magnetization in the z-direction. The ground state corresponds to the ferromagnetic spin 
state 



1 N I \ 1 

M = ^EKi> + < s 2 Vj = 5 
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Altogether, it could be concluded that the IF anisotropy in the Heisenberg intra-rung coupling is re- 
sponsible for quantum features of otherwise classical SR and FM states and, moreover, it may also lead 
to the appearance of a new disordered phase QPM2. Two paramagnetic phases QPM1 and QPM2 have 
quite similar features: they are both disordered states with the energy gap in their excitation spectrum 
and, consequently, their pair spin-spin correlation functions decay exponentially. Both quantum param- 
agnetic phases can be distinguished by the square of zth component of the total spin (S?) 2 on iih rung, 
which is equal to (Sp 2 = 0(1) in QPM1 (QPM2). 

Now, let us pay our attention to the ground-state phase diagram established for the particular case of 
XZ-Ising ladder as depicted in figure [2]by considering ]\ = 0. Assuming the XZ intra-rung interaction, 
one gets a striking competition between the X-X interaction along the rungs and the Z-Z interaction 
along the legs, while the additional Z-Z interaction along rungs acts generally against the X-X inter- 
action. It should be also mentioned that one may recover some known examples from the ground-state 
phase diagram of the XZ-Ising ladder presented in figure[2] In fact, figure[2](a) shows the particular lim- 
iting case of a quantum compass ladder with an additional diagonal frustrating Ising interaction. Let us 
follow the known results of a simple quantum compass ladder [3] to be obtained from our model by dis- 
regarding the frustrating Ising interaction J% - 0. Both Hamiltonians H + and H~ (23) become identical 
under this special condition and, consequently, the ground state of the model is always two-fold degen- 
erate due to the equality = . The investigated model is either in the SL or FM state for Ijf < -Jz, 
either in QPM1 or QPM2 state for \ jf> \h\, either in SR or Neel state for \ jf< h- Note that the quantum 
phase transition from the disordered to the long-range ordered state takes place at \ Ji | = \ jf . 

It is quite evident that the diagonal interaction h removes the two-fold degeneracy of the ground state 
[see figure[2](a)] by changing the effective spin interaction in the effective Hamiltonians H + and H~ . The 
effect of the additional Z-Z intra-rung interaction is shown in figure [2] (b)-(d), where the results for a 
different relative strength of Jf are presented. It could be understood from <2.5> that Jf lowers initially 
the energy of the . = -\ configuration by the amount \NJ* with respect to the energy of the . = 
| configuration. Therefore, the regions of the phases, which correspond to the 5j ,• = | configuration, 
become smaller with increasing Jf. For sufficiently strong Jf, QPM2 phase can completely disappear 
from the ground-state phase diagram. Another distinctive feature is that Jf makes the singlet-dimer phase 
favorable along the line Jz - h inside the region of QPM1 state. 

The relevant ground-state behavior can be supported by the dependencies of respective order pa- 
rameters as displayed in figure [3] (a). If there is no frustrating diagonal interaction, the model may stay 
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Figure 2. (Color online) Ground-state phase diagram of the XZ-Ising ladder for jf = 1, Jz = and four 
different values of /J: (a) ]\ = 0; (b) /* = 0.05; (c) Jf = 0.1; (d) Jf = 0.5. QPMl is the paramagnetic phase 
with the prevailing rung states \cj)' Q Q ). QPM2 is the paramagnetic phase with the prevailing rung states 
\(p\ _>. Dotted lines indicate the rung singlet-dimer state. 
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in the disordered QPM1 or the ordered SL and Neel phases. The corresponding nearest-neighbor corre- 
lation function along the legs shows a continuous change with a weak singularity at the quantum critical 
points indicated by filled circles in figure[3](b). The curve for another particular case J$ - 0.25 looks sim- 
ilar except that the diagonal interaction of this strength leads to a direct phase transition between two 
disordered quantum paramagnetic states QPM1 and QPM2. This unusual transition can be recognized 
from the relevant dependence of the nearest-neighbor correlation function, which sustains a jump at this 
special critical point. It is interesting to note that the further increase of a frustrating Ising interaction 
Js demolishes both disordered phases QPM1 and QPM2. Thus, one may also detect the quantum phase 
transition between two ordered SL and Neel phases, whereas the order parameters do not reach zero 
continuously in this particular case. 




In flgure|U the ground-state phase diagram of the XF-Ising ladder is depicted by considering another 
particular case with Jf - 0. The effect of the Y-Y intra-rung interaction has some similarities with the one 
of Z-Z intra-rung interaction, although the origin is completely different. The ground-state energy of H~ 
is generally lowered with respect to that of H + , because ft > increases the effective transverse field in 
H~ and decreases it in H + . Owing to this fact, the regions of SL, QPM1, and Neel states that correspond to 
the uniform configuration with all 5^ . — — | get extended with an increase of ]\. One may also generally 
conclude that ] x > acts against the QPM2 phase, which gradually shrinks with increasing /J until it 
completely disappears at /j = Similarly to the case with the Z-Z intra-rung interaction, the Y-Y 
intra-rung interaction also induces the presence of the rung singlet-dimer state along a special line — h- 
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4. Conclusions 

In the present paper, the effect of the most general XYZ anisotropy in the intra-rung interaction 
on the ground state of the spin-i Heisenberg-Ising two-leg ladder was investigated in detail. It has been 
shown that the most general kind of anisotropy, which breaks the rotational symmetry of the Heisenberg 
interaction, may lead to the appearance of new quantum phases in the ground-state phase diagram. We 
have also considered the special case of quantum compass ladder with an additional frustrated diagonal 
interaction and showed that the singlet-dimer phase cannot appear in this particular case. The order 
parameters and the nearest-neighbor correlation function were calculated and analyzed in detail in the 
ground state. It has been demonstrated that the relevant behavior of the correlation function can help us 
to reveal the quantum phase transition between two different disordered quantum paramagnetic phases. 
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OCHOBHMM CTdH cniH-1/2 flBOHOrOl' flpa6l/IHKI/l 

raM3eH6epi i a-l3i/iHi'a 3 XYZ B3aeM0flieK> B3flOB>i< ma6/iiB 

T. Bepxcviflhp] M. CrpeHKaP 

1 iHCTMTyT (J>i3HKl<l KOHfleHCOBaHWX CHCTeM HAH YKpaiHH, By/1. CBEHUjUbKOrO, 1, 7901 1 ZlbBiB, YKpaiHa 

2 Kacf>e,qpa TeopeTUHHOi (f>i3WKH i anpo(|)i3HKH, YHiBepcuTeT n.M. LUaefiapwKa, napK AHre/imyM, 9, 04001 Kouinuj, 
OiOBaubKa pecny6/iiKa 

KBamoBy cniH-1/2 flBOHOry flpa6wHKy 3 aHi30TponHOK) XYZ B3aeMOflieio raw3eH6epfa B3flOBX ma6/iiB i I3HH- 
r'OBOKj B3aeMOflieK) Mix cniHaMH Ha cyciflHix ma6/iax po3r/iHHyTO b Mexax CTpororo niflxo/ry, hkum fpyHTyeTbcn 
Ha yHiTapHOMy nepeTBopeHHi. MacTKOBnw BwnaflOK uopfim 3 X-X B3aeMOflieio b3,qobx Lua6/iiB Bi,qo6paxae 
MOfle/ib KBaHTOBOro KOMnacy Ha flpa6wHuj 3 flOflaTKOBWMn fliaroHa/ibHUMW IsnHroBWMn B3aeM0flinMn. Buko- 
Phctobk)4w yHiTapHe nepeTBOpeHHH, MOfle/ib MOXHa 3Bec™ flo nonepenHoro /iaHU_K)XKa l3HHra 3 KOMno3HT- 
hmmh cniHawiw, i ak HaaiiflOK ochobhuw CTaH MOXHa 3HaPuw CTporo. Mn OTpnMyeMO c|>a30By fliarpaiviy ochob- 
Horo CTaHy i aHa/ii3yeMO B3a€MOBn/inB KOHKypeHuji Ki/ibKOx cfiaKTopiB: XYZ aHi30Tponii B3aeMOflii Tati3eH6ep- 
fa, I3hhi , oboi B3aeMOflii B3flOBX ziaHuroxKiB Ta cfipycTpoBaHOi fliaroHa/ibHOi ISHHroBOi B3aeMOflii. Po3i7iflHyTa 
MOfle/ib flewioHCTpye CK/iaflHi (f>a30Bi fliarpaMW ochobhoto CTaHy, BK^KD4aKJ4n fleKi/ibKa HemnoBnx KBaHTOBnx 
BnopaflKOBaHnx CTaHiB, flBa pi3Hwx HeBnopaflKOBaHnx KBaHTOBnx CTaHiB, a TaKox KBaHTOBi cf>a30Bi nepexoflu 
nepuioro i flpyroro nopaflKy Mix 14MMH cf>a3aMW. 

K/noMOBi c/ioBa: mamoBa cniH-1/2 Apa6nHKa, tohhi pe3ynbTam, ochobhi/iH craw, KBaHTOBMti tpa30Bnii 
nepexifl 
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